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ABSTRACT 

Spin parameters of stellar-mass black holes in X-ray binaries are currently 
being estimated by fitting the X-ray continuum spectra of their accretion disk 
emission. For this method, it is necessary to know the inclination of the X-ray- 
producing inner region of the disk. Since the inner disk is expected to be oriented 
perpendicular to the spin axis of the hole, the usual practice is to assume that 
the black hole spin is aligned with the orbital angular momentum vector of the 
binary, and to estimate the inclination of the latter from ellipsoidal modulations 
in the light curve of the secondary star. We show that the inclination of the disk 
can be inferred directly if we have both spectral and polarization information 
on the disk radiation. The predicted degree of polarization varies from 0% to 
5% as the disk inclination changes from face-on to edge-on. With current X-ray 
polarimetric techniques the polarization degree of a typical bright X-ray binary 
could be measured to an accuracy of 0.1% by observing the source for about 
10 days. Such a measurement would constrain the disk inclination to within a 
degree or two and would significantly improve the reliability of black hole spin 
estimates. In addition, it would provide new information on the tilt between 
the black hole spin axis and the orbital rotation axis of the binary, which would 
constrain any velocity kicks experienced by stellar-mass black holes during their 
formation. 

Subject headings: accretion, accretion disks — black hole physics — polarization 
— radiative transfer — relativity — X-rays: binaries 

1. Introduction 

X-ray polarimetry is an uncharted frontier with great promise for the exploration of the 
behavior of accreting black holes and for the measurement of black hole spin. To date, only 
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small X-ray polarimeters have been flown and the only result of consequence has been the 
measurement of the polarization of the Crab Nebula (Novick et al. 1972; Weisskopf et al. 
1978). Anticipating the launch of larger and more sensitive X-ray polarimeters in the future, 
this paper discusses possible uses of polarization observations to study thermal emission from 
accretion disks around stellar-mass black holes in X-ray binaries. 

Polarization measurements at radio and optical wavelengths attest to the power of 
polarimetry. In radio astronomy, the measurement of polarization is trivial. However, its 
consequences have been of central importance since Alfven & Hcrlofson (1950) first suggested 
the synchrotron mechanism as the source of celestial radio emission. For example, studies 
of this strongly polarized radiation has allowed the mapping of magnetic fields that are 
threaded throughout stellar coronae, supernova shocks, relativistic jets, and the interstellar- 
and intergalactic media. At optical wavelengths, the measurement of polarization is more 
difficult, but the payoffs have also been large. For example, polarimetric observations of 
NGC 1068 provided a spectacular confirmation of the unified model of AGN (Antonucci 
& Miller 1985). Polarization measurements are most challenging at X-ray wavelengths; 
however, sure instrumental approaches are known (§6), and the discovery space and potential 
rewards are large. 

X-ray polarimetric information - direction and degree - increases the parameter space 
used to investigate compact objects from the current two - spectra and time variability - 
to four independent parameters that models must satisfy (Rees 1975; Lightman & Shapiro 
1975; Mcszaros et al. 1988). X-ray polarimetry can provide qualitatively new information 
on source geometry and magnetic fields on spatial scales comparable to a black hole event 
horizon. In this paper, we follow the pioneering work of Stark & Connors (1977), Connors 
& Stark (1977), and Connors et al. (1980), and we narrow our focus to consider solely the 
problem of polarized thermal emission generated by a thin accretion disk around a stellar- 
mass black hole. The theory developed here of course applies as well to thermal accretion 
disks around supermassive black holes. We do not consider coronal effects that can generate 
an X-ray refiection spectrum (e.g., Ross et al. 1999; Dovciak et al. 2004) and Comptonized 
components of emission (e.g., Gierlihski et al. 1999). 

The polarization model presented here extends our workhorse accretion disk model 
KERRBB (Li et al. 2005), which we have used to estimate the spins of four stellar-mass 
black holes via fits to their thermal continuum spectra (Shafee et al. 2006; McClintock et al. 
2006; Liu et al. 2008). The extension of KERRBB presented herein includes a simple electron 
scattering atmosphere based on an a- viscosity prescription of the stress (Shakura & Sunyaev 
1973) and a treatment of the polarized radiation field via the Stokes parameters. Radiative 
transfer effects will be considered in a future paper. 
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In §2, we discuss our motivations for this work, focusing in particular on the apphcation 
of polarization measurements for estimating the spins of accreting black holes. In §3, we 
briefly review the polarization of disk emission induced by electron scattering in an accretion 
disk atmosphere, assuming a semi-infinite plane-parallel medium. In §4, wc present the 
mathematical scheme for calculating the polarization as observed by a remote observer, 
after the polarized disk emission propagates through the curved spacetime of a Kerr black 
hole. In §5, we apply our polarization calculations to the measurement of black hole spin, 
and in §6, we briefly discuss the prospects for observing the predicted polarization signal 
from accreting stellar-mass black holes. Finally, in §7, we summarize the results. Some 
mathematical details related to the thin disk model we employ and the propogation of the 
polarization vector in Kerr spacetime are presented in Appendices A and B. 



2. Motivations 

There are two principal and related motivations for undertaking this work: (1) to val- 
idate a key assumption of the continuum-fitting method of determining spin, and (2) to 
determine if the spins of stellar mass black holes are aligned with the orbital plane. We 
discuss these topics in turn in the following subsections. We then conclude this section by 
commenting on three additional avenues for determining black hole spin, viz., the Fe K line, 
high-frequency QPOs and X-ray polarimetry. 



2.1. Validation of the Continuum-Fitting Method of Meeisuring Spin 

Using the continuum-fitting method, we have published spin estimates for four stellar 
black holes (§1), and wc expect to publish spin estimates for several more during the next 
few years. In this work, it is essential to use thermal dominant (formerly high soft state) 
data (McClintock & Rcmillard 2006). A feature of this state is the dominance of a soft 
blackbody-like component that is emitted by optically-thick gas in the accretion disk. A 
minor nonthermal tail component of emission is also often present, but it contributes typically 
only a few percent of the ~ 1 — 10 keV flux. Thus, these thermal-state spectra are largely 
free of the uncertain effects of Comptonization and are beheved to match closely the classic 
thin accretion disk models of the early 1970s (Shakura & Sunyaev 1973; Novikov & Thorne 
1973). There is an abundance of such suitable thermal-state continuum spectra available in 
the NASA HEASARC archives for most black hole binaries. These spectra are ideally suited 
for study using the polarization model presented in this paper. 
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The promise of the continuum-fitting method for determining black hole spin is attested 
to by 20 years of experience that has demonstrated for the thermal dominant state the 
presence of a stable inner disk radius (see §4 in McClintock ct al. 2008). The inner edge is 
presumably closely tied to the innermost stable circular orbit (ISCO). X-ray observations, 
coupled with ground-based optical and infrared observations, are used to measure the radius 
of this orbit. Applying the continuum-fitting method is analogous to the familiar problem of 
determining the radius of a star of known distance given its effective temperature and flux 
(McClintock et al. 2008). Accordingly, this method requires knowledge of the luminosity of 
the X-ray source, which depends not only on its distance but also on the projected area and 
hence inclination of the accretion disk. Accurate values of the inclination of the orbital plane 
iorb, distance D, and black hole mass M are readily obtained from ground-based observations. 
This optical work has been underway for 35 years and is now undergoing a renaissance (e.g., 
Orosz et al. 2007). 

The continuum-fitting method is straightforward to apply. All the required data are 
readily obtainable, and even the theory of disk accretion in strong gravity is tractable 
(Narayan et al. 2008; Shafee et al. 2008b). However, the reliability of the continuum- 
fitting method is called into question by a single assumption, viz., that the plane of the inner 
X-ray-emitting portion of the disk is aligned with the binary orbital plane, whose inclina- 
tion angle iorb is determined from optical observations. Unfortunately, the continuum-fitting 
method cannot fit for the inclination of the inner disk and check for disk warp (§2.2) because 
there is a degeneracy between the inclination and spin parameters (§ 5.2; see also Liu et 
al. 2008). As we show, polarimetry observations can easily determine the inchnation of the 
inner disk. This is our principal motivation for writing this paper. 



2.2. The Inclination of the Inner Accretion Disk 

Whether a black hole's spin is aligned with its accretion disk and the time scale for 
any misalignment to dissipate are topics that have received considerable attention recently 
(Fragile & Anninos 2005; King et al. 2005; Lodato & Pringle 2006; Fragile et al. 2007; Martin 
et al. 2007). In this paper, we restrict our attention to the question of alignment in the case 
of stellar black holes in X-ray binaries. If a black hole's spin were to be misaligned from the 
orbital vector, the inner disk would be warped away from the outer disk by Lense-Thirring 
precession and the Bardccn-Pcttcrson effect (Bardccn & Pcttcrson 1975), which would be 
expected to cause a global precession of the main body of the disk and have important 
observational consequences (e.g., enhanced mass accretion rate and jet precession; Fragile et 
al. 2007; Martin et al. 2007). 
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Are significant misalignments expected in the case of accreting stellar black holes? On 
the one hand, it seems plausible that the spin will be aligned with the orbit vector given the 
strong tidal forces acting on the pre-supernova star (Zahn 1977) and, as generally assumed 
(e.g.. Fryer & Kalogera 2001), that collapse of the stellar core to a black hole will be nearly 
spherical and mostly free of the strong shocks and velocity kicks that nascent neutron stars 
suffer. Even if shocks and kicks do occur, it is expected that their effects will be less for black 
holes by the factor Mbh/A^ns ~ 7 (e.g.. Fryer & Young 2007). On the other hand, there 
is some evidence for warped disks based on radio-jet data (Maccarone 2002), although this 
evidence is weak (Narayan & McClintock 2005). One can attempt to infer inchnation via 
studies of radio jets (Mirabel & Rodriguez 1999) and Fe K emission line profiles (Fabian et 
al. 2000). However, as we show in this paper, the direct and sure approach is by measuring 
the degree of polarization of the thermal accretion-disk radiation, which depends strongly 
on the inclination angle idisk of the inner disk. 



2.3. Three Additional Methods of Meeisuring Black Hole Spin 

The Fe K-line method of determining spin is applicable to supermassive black holes as 
well as stellar black holes (Miller 2007; Reynolds & Fabian 2008). It has received considerable 
attention and has recently yielded quantitative estimates of spin for the Seyfert 1.2 galaxy 
MCG-6-30-15 (Brenneman & Reynolds 2006) and for the steUar black hole GX 339-4 (Miller 
et al. 2008). 

High-frequency QPOs may develop into a precise method of determining spin once the 
correct model of the oscillations has been identified (e.g.. Wagoner et al. 2001; Torok et 
al. 2005; Schnittman 2005). X-ray polarimetry may also provide a means of determining 
spin. This possibility is implicit in the work of Connors et al. (1980), which shows that for 
reasonable assumptions one expects the plane of polarization to swing smoothly through a 
large range of angle as the photon energy increases. By modehng this effect and the degree 
of polarization one can hope to obtain an independent estimate of black hole spin, a topic 
that we touch upon in §5.1. 

Thus, there are potentially a total of four methods of measuring black hole spin, two of 
which - the continuum-fitting method and the Fe K-line method - are presently delivering 
results. Because black hole spin is such a fundamental parameter, it is important to attempt 
to measure it by as many of these methods as possible, as this will provide arguably the best 
possible check on the results. 
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3. Polcirization of Disk Emission Induced by Electron Scattering 

We discuss in this and the following section the technical details of how we calculate the 
polarization of the thermal disk emission. We present additional details in the Appendices. 
The presentation is pedagogical so that this paper and Li et al. (2005) together provide a 
complete and self-contained description of the calculations. 

The radiation generated deep inside an accretion disk is initially unpolarized, but it 
becomes polarized as a result of electron scattering in the disk atmosphere. For an optically 
thick and geometrically thin accretion disk, the disk atmosphere is locally well described by 
a semi-infinite plane-parallel medium. We make this approximation. 

A radiation field is completely described by four Stokes parameters: /(=// + /,., the 
total intensity), Q {= Ii — Ir), U, and V (Chandrasekhar 1960). Here the subscripts / and 
r refer, respectively, to polarization in the meridian plane (i.e., the plane containing the 
symmetry axis and the line-of-sight to the observer) and at right angles to it. In a plane- 
parallel atmosphere with no incident radiation, the axial symmetry of the radiation field 
requires that the plane of polarization be either along I or r. Therefore, U — V — 0, and the 
two parameters / and Q, or equivalently the two intensities Ii and Ir, suffice to characterize 
the radiation. The degree of polarization of the radiation field is then given by 

p^\{Q^ + u' + vY'^^-^;^. (1) 

At the surface of a pure scattering atmosphere, the two intensities // and are deter- 
mined by 

hil,) = -^FHi{ix)q , Irili) = -^FHr{li){li + c) , (2) 

where /i = cos^, 9 is the polar angle of the direction of propagation of the photon to the 
disk normal, F is the total emitted flux density. Hi and are Chandrasekhar's //-functions 
(Chandrasekhar 1960), and q and c are two integral constants. Following Chandrasekhar 
(1960), we have suppressed the dependences of the relations on the photon energy since no 
ambiguity is fikely to arise. The constants q and c are given by (Chandrasekhar 1960, §X) 

8(A + 2ai)-6(Aoai + aoA) A q^\"^ 

^ = 2,{Al + 2al) ' '^V-^) ' 

where q;„ and are the moments of order n of Hi{^) and Hr{ii): 

an= [ Hi{iJi)iJ'diJi , Ar,= I Hr{iJi)iJ'diJi (n > 0) . (4) 
JO Jo 
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Figure 1 shows the degree of polarization P as a function of ji for a scattering atmo- 
sphere. At /i = 1 = 0), we have P = because of symmetry. As ^ decreases (i.e., 9 
increases), P increases monotonically and reaches its maximum value of 11.7% at = 
{9 = 7r/2). For a semi-infinite atmosphere we always have Ir{^) > (sec table XXIV 

in Chandrasekhar 1960), except at ;U = 1 = 0) where /r(/u) = /«(0,/i). Hence we have 
ip = t:/2 for all /x, where ip is the angle of the linear polarization vector relative to the 
meridian plane. 

The above results arc for a pure scattering semi-infinite atmosphere. However, an accre- 
tion disk also has absorption processes which may be important, especially at large disk radii 
where the temperature is low. Absorption tends to destroy the polarization of photons and 
hence to reduce the degree of polarization. We defer a detailed computation of the effect of 
absorption to a later paper. Here, following Laor et al. (1990) and Chen & Eardley (1991), 
we simply reduce the degree of polarization by the following approximate factor 

?w = — ^ — , (5) 

where is the optical depth due to electron scattering, and Tab is the optical depth due to 
photon absorption. Thus, we write the degree of polarization of photons emerging from the 
disk surface as 

-Pem — QwPerafl , (6) 

where Pcm,o is the degree of polarization for a pure electron scattering atmosphere and is 
given by equations (1) and (2). We assume that the orientation of the polarization is not 
affected by absorption. 

To complete the description of the polarized emission from an accretion disk, we need 
the flux F and the optical depths t^s and Tab as functions of the radius R. These are obtained 
using the disk model described in Appendix A. 



4. Poleirization of the Disk Emission as Observed by a Remote Observer 

An arbitrary radiation field can be decomposed into an unpolarized component, with 
Stokes parameters 

{(1 - P)7, 0, 0, 0} , (7) 
and a completely polarized component, with Stokes parameters 

{PI, Q, U, V} , (8) 
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where P is the degree of polarization. Further we have 

where ip is the angle between the major axis of the polarization ellipse and the direction 

described by the sTibscript / (see §3)/ and tan/3 is the ratio of the minor axis to the major 
axis of the polarization ellipse. Photons emitted by the disk are linearly polarized, so we 
have V — and /3 = 0. Then, by the definition of the degree of polarization P, we have 

g = /pC0s2V', C/ = /pSin2V', Ip = PI- (10) 

Equivalently, we can write 

Q + tU = Ipe^'"^ . (11) 

When several independent streams of hght are combined, the Stokes parameters for 
the combined radiation are the sums of the respective Stokes parameters of the individual 
streams (Chandrasekhar 1960); we refer to this as the superposition theorem. Therefore, we 
may consider the unpolarized and polarized components of the radiation separately. 

For the polarized component, consider a beam of perfectly polarized radiation emitted by 
an infinitesimal surface element on the disk and received by an observer along an infinitesimal 
solid angle element dflobs- Let this radiation have intensity /p,obs and polarization angle ipobs 
as measured by the observer. Then, the observed Stokes parameters are 

{-^p,obs, Qobs, t^bs, 0} , (12) 

where 

Qohs = Ip,ohs cos 2'0obs , ^obs = -^p.obs sin 2V^obs ■ (13) 

Summing over the radiation received from all disk elements we have, by the superposition 
theorem, 

(-^p.obs) = / -^p,obsC?^obs , (14) 



obs 



and 



(Qohs) + i{Uohs) = -To / (Qohs + iUohs) d^ohs 

^"obs J 

= -TT^ / /p,obse'*'^°''=ciOobs , (15) 

^"obs J 



^We denote the angle of the plane of polarization from the direction I by ip (rather than the % in 
Chandrasekhar 1960) to be consistent with Connors et al. (1980). 
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where '( )' denotes an average over solid angle, and AQobs is the total solid angle subtended 
by the disk on the sky. 

Note that, in general, we have 

(Qobs)'+(t/obs)'<(/p,obs)', (16) 

i.e., the radiation is not fully polarized. Even though the individual contributions from each 
disk element may be perfectly polarized, the radiation received by the observer can still 
become partially polarized through averaging because the individual rays may have different 
values of ■0obs at the observer. In the context of an accretion disk around a black hole, even if 
all the emitted rays have the same polarization direction at their points of emission, geodesic 
propagation can cause changes in i/j along each photon trajectory and cause a reduction in 
the observed degree of polarization. In other words, the rotation of the disk and the spin 
of the black hole can alter the polarization state of the original radiation and destroy the 
observed polarization at some level. (In contrast, whatever unpolarized radiation is emitted 
by the disk remains unpolarized at the observer; thus, geodesic propagation and averaging 
can only decreased the degree of polarization.) 

For the unpolarized component, we have 

(-^u,obs) = -Tri / -^u,obsC?^obs , (17) 



A^obs _ 

where /u,obs = -^obs ~ -^p.obs- By equations (14) and (17) we have 

(-^u.obs) + (-^p.obs) — -To / (-^u,obs + -^p,obs) O^^^obs = Xri / -^obsO^^obs = (-^obs) , (18) 

^"obs J ^"obs J 

i.e., the total intensity is conserved. 

Using the fact that Iei^^/E^^^ is invariant along the path of a photon, where £"100 is the 
photon energy measured by a local observer along the photon path (Misner et al. 1973), 
equation (15) can be rewritten as (c.f. eq. 9 of Connors et al. 1980) 

(Qobs) + ^(^obs) = I ^'Pem/eme'*^°'^=(ifiobs , (19) 

^"obs J 

where Pem is the degree of polarization of the radiation at the time that the radiation emerges 
from the disk surface, and g is the redshift factor of the photon (Li et al. 2005). Similarly, 
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we can rewrite the total intensity as' 



2 



(/obs) = -TTS— / g'leradn^hs ■ (20) 



AQobs 

The observed average degree of polarization is then given by 



{P)^j^V(Q^^^f + {U^ , (21) 

\-'obs/ 



3 



while the observed average angle of polarization, (ip), is determined by 

sin(2(V')) = , = , cos(2(V')) = , = . (22) 

v/(Qobs)2 + (f/obs)2 V(Qobs)^ + (f/obs)2 

In the definitions of the above dimensionless quantities the total solid angle AQobs cancels 
out. 

The solution to equation (22) is 

(^) = (^)pr + n7r, (23) 
where n is any integer, and the primitive angle (■0)pr (defined to be in the range O-tt) is given 

by 

/./A ^! 1 arccos^Q , ((t^obs) > 0; < (^)pr < f ) , . 

^'f'hr I ^ _ 1 ^^^^^3^^ ^ (^f^^^^^ ^ Q. ^ ^ ^^^^^ ^ , I J 

where 

^ ^ . (25) 

V(Qobs)2 + (t/obs)2 

In equation (19), the degree of polarization Pem is evaluated at the disk surface by equation 
(6), but the angle of polarization ipohs is evaluated at the observer and is related to the angle 
at the disk surface by the propagation equation of the polarization vector. The propagation 
equation is described in detail in Appendix B. 



"^Note that, in order to simplify the notation, we have suppressed the dependenc on £^em - the photon 
energy as measured by a local observer corotating with the disk - in our expressions. This should cause no 
ambiguity. 

^We do not use tan(2(i/))) = (C/obs)/(Qobs)- This single equation cannot determine the value of (tp), since 
tan(2(V')) is unchanged under the transformation (Qobs) ~* ~(Qobs) and (f/obs) ~* —{Uobs)- This is related 
to the fact that the primitive period of sin a; (and cos a;) is 27r while the primitive period of tana; is tt. 
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In the limit of a semi-infinite plane atmosphere, the polarization vector of the emerging 

radiation is in the disk plane, so wc have ipcm = 7r/2 (see §3). Since V'obs appears in equation 
(19) in the form of cos2'?/'obs and sin2^obs) addition of rni to tpobs (^ — =tl) ±2, ...) does not 
change the final results. Hence, at the observer, by equation (B74) we can write 

V'obs = - ^GR , (26) 

where $gr is calculated by equations (B69), (B70), and (B75) as described in Appendix B. 

Relativistic effects on photon polarization are reflected in the above formula in two 
apsects. First, relativistic effects lead to a rotation to the plane of polarization by an angle 
$GR- Second, because of hght bending the angle of the photon wave vector from the disk 
normal as the photon leaves the disk surface, 9, differs from the disk inclination angle, idisk- 
The Newtionian limit is obtained by setting $gr — and 6 = i^isk- By equation (26) we 
then have ■0obs — 7i'/2. Then, by equation (19) we have (L^obs) — and 

(Qobs) = - f^T'" / qwg^Ierad^ohs , (27) 
^"obs J 

where equation (6) has been used, and Pem,o — Pem,o{& — ^disk)- Hence, in the Newtonian 
limit we have (■0) = 7r/2, and 

The high energy spectrum is dominated by photons emitted by the inner disk region where 
photon absorption is not important {q^ ~ 1), hence we have (P) ~ -Pem,o for high energy 
photons. 

As we will see in the next section, relativistic effects make the observed features of 
photons emitted by the disk dramatically different from that predicted by the Newtonian 
theory. 



5. Applications to the Measurement of Black Hole Spin 

5.1. Compcirison with Connors, Piran Steirk (1980) 

Since the present work is a modern update on the pioneering work of Connors et al. 
(1980), we begin by recomputing the models discussed in their paper. Appendix B describes 
in detail how we calculate the polarization angle of an inflnitesimal photon stream as observed 
by an observer at inflnity. A key quantity is the primitive rotation angle $gr deflned by 
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equations (B69), (B70), and (B75). The calculation of $gr requires an evaluation of the 
dimensionless quantities X, Y, S, and T (eqs. B79-B81), in addition to the orientation of 
the photon wave vector relative to the disk plane {6 and 0), the photon redshift factor g, 
and the Lorentz factor F of the disk fluid. All these quantities can be evaluated with the 
formulae given in Appendix B of this paper and Appendix C of Li et al. (2005). 

The polarization angle i/^obs is related to the rotation angle $gr by equation (26). After 
calculating the angle ■?/'obs for each light ray that reaches the observer within an infinitesimal 
solid angle dQ^hs, one can evaluate the integrals in equations (19) and (20) with the ray- 
tracing technique described in Li et al. (2005). Then the net degree of polarization (P) 
and angle of polarization (ip) of the combined radiation, as seen by the observer, can be 
calculated via equations (21) and (23)-(25). In the calculations reported here, the effect 
of photon absorption has been taken into account with the simple approach outlined in §3. 
That is, photon absorption reduces the degree of polarization by introducing a multiplying 
factor q^f defined in equation (5). 

Figure 2 shows the variation of the observed polarization angle {ip) and the observed 
degree of polarization (P) as a function of photon energy for two disk inclinations idisk and 
various values of the dimensionless spin parameter a* = cJ/GM^, where J is the angular 
momentum of the black hole. All the models shown in Fig. 8 of Connors et al. (1980) 
are included here; several other models are also shown for completeness. Qualitatively, our 
results arc similar to those obtained by Connors et al. (1980). Quantitatively, however, 
there are unexpectedly large differences. Compared to Fig. 8 of Connors et al. (1980), our 
results in Fig. 2 predict less photon absorption at low photon energies (indeed this can be 
seen clearly in the Newtonian hmit, the dotted lines in the figure). The results for high 
photon energy are also different.^ 

Rather surprisingly, the models differ considerably even in the basic shape of the contin- 
uum spectrum, as seen by comparing Fig. 7 of Connors et al. (1980) with the corresponding 
results from our model shown in Fig. 3. We have assumed that the disk emission is purely 
blackbody, and we have allowed for the spectral hardening due to electron scattering and 
Comptonization through a constant hardening factor f^oi- The spectra shown by Connors 
et al. do not appear to be pure blackbody (see their Fig. 7). Indeed, their model seems 



^Note that we have set M = 7 x lO^^gs"^ in the models shown in Figs. 2 and 3, since this is the value 

mentioned by Connors et al. (1980) in the second sentence of § IVc. In their caption to Fig. 5, however, 
they mention M = 10^''gs~^. We have computed models with the latter value of M (not shown) and the 
results are nearly the same. 
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to have too much emission both at low and high energies.^ Since we do not have sufficient 
information on the precise assumptions made by Connors et al. (1980) in their calculations, 
we are unable to resolve the discrepancy. 

We have also compared our results with the calculations reported by Agol & Krolik 
(2000) and Dovciak et al. (2008). The easiest and most direct comparison is when returning 

radiation and photon absorption are neglected, which has been assumed in the model of 
Dovciak et al. (2008), amd corresponds to the dashed lines in Fig. 11 of Agol & Krolik 
(2000). For this case, we find very good agreement with the results from our code for all 
quantities of interest: flux, degree of polarization, and polarization angle. When returning 
radiation is included, Agol & Krolik (2000) calculute the polarization induced by scattering 
off the surface of the disk. This effect is currently neglected in our model. We also neglect 
the effect of Faraday rotation due to magnetic fields in the disk atmosphere (Agol et al. 
1998). 

Before concluding this subsection, we would like to note the following interesting point. 
Figure 2 shows that, for a fixed M and idisk, the polarization of the observed radiation 
exhibits a large variation as a function of the spin parameter. Therefore, if we have already 
measured M and i^isk (see the next two subsections), then we could use polarization data 
to constrain a*. This would provide an estimate of a* independent of the disk continuum 
fitting method. In fact, with sufficiently high quality polarization data we might even be 
able to solve for both ijisk and a* without using the continuum spectrum at all. 

5.2. Spin-Inclination Degeneracy in the Disk Continuum Spectrum 

We proceed now to discuss the continuum fitting method and the important role that 
polarization measurements could play in this method. To determine the spin parameter a* 
of a black hole in an X-ray binary, we fit the X-ray continuum spectrum of its accretion 
disk using the publicly- available model KERRBB. Three parameters of the system need to 
be measured independently: the mass of the black hole M, the distance to the black hole 
D, and the disk inclination angle idisk- kerrbb has two principal fit parameters: the disk 



^At low photon energies, the blackbody photon flux density should have the well-known asymptotic form 
Ne oc E~'^/^, e.g., our Fig. 3, or equivalently the energy spectrum should go as Fe oc E^/^. However, Fig. 7 
of Connors et al. (1980) seems to behave as Ne oc E^'^ near E = 0.1 kcV. Also, their model C, which 
corresponds to a rapidly spinning black hole (a, = 0.998) has an imphysically hard spectrum with emission 
extending up to 100 keV, which is much too hot for thermal emission from a thin accretion disk around a 
stellar-mass black hole. 
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mass accretion rate M and the dimensionless spin parameter a*. The mass accretion rate 
serves the role of a normahzation factor and depends on the X-ray luminosity, while a* is 
determined principally by the temperature or hardness of the X-ray spectrum. Generally, 
it is straightforward to estimate these two parameters from X-ray continuum data. In fact, 
even when the spectral data require additional features for a good fit, e.g., photoelectric 
absorption, an iron line, a weak Comptonized tail, one can still estimate M and a* with 
good precision. 

In addition to M, D and idisk, the spin determination requires a theoretical model of 
the disk flux profllc F{R) as a function of radius R, and a detailed disk atmosphere model to 
estimate the spectral hardening factor ff.o\- For F{R), kerrbb makes use of the relativistic 
disk model of Novikov & Thorne (1973). The no-torque inner boundary condition assumed 
in this model has been questioned (e.g., Krolik 1999; Krolik & Hawley 2002; Krolik et al. 
2005), but recent work suggests that it may not be a serious issue when the disk luminosity 
is low and the disk is geometrically thin (Paczyhski 2000; Li 2002; Afshordi & Paczyhski 
2003; Shafee et al. 2008a,b). With regard to /coi, disk atmosphere models including metal 
opacities have been recently computed and /coi has been estimated as a function of the disk 
luminosity and inclination (Davis et al. 2005; Davis & Hubcny 2006; Davis et al. 2006). 
This is a significant improvement over previous models which included only free-free opacity 
(Shimura k Takahara 1995). 

We thus expect that a reliable determination of a* of a stellar black hole can be made 
when we have high quality X-ray data in the thermal state — usually not a problem — when 
we have good measurements of M, D and idisk- For a number of stellar black holes M and 
D are well measured (Charles & Coe 2006; Orosz 2003; Orosz et al. 2007). However, the 
inclination angle of the inner disk cannot be directly measured, except in a few cases when 
the system has a radio jet and the orientation of the jet (which is presumably perpendicular 
to the inner disk plane) can be determined (Hjellming & Rupen 1995; Hjellming et al. 2000; 
Orosz et al. 2001). For the remaining systems, what we can measure is the inclination of 
the binary orbit iorb, and we have to assume that the inner disk is aligned with the binary 
orbital plane: idisk = "^orb- This assumption is probably reasonable (see § 2.2), but it would 
be better if we could avoid it altogether. 

In view of the uncertainty in the disk inclination, could we obtain idisk directly from 
X-ray spectral data? Unfortunately, this is not possible since a given observed continuum 
spectrum can be fitted well with different combinations of and idisk- Table 1 shows an 
example of a set of degenerate models with distinct combinations of a* and idisk- The mass of 
the black hole is fixed at lOM©, the distance to the black hole is fixed at 10 kpc, the spectral 
hardening factor is assumed to be fixed at 1.6, and the mass accretion rate M of each model 
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has been adjusted such that all the models have the same luminosity. Figure 4 shows the 
disk continuum spectra of these models. Although the models have very different values of 
and idisk, we see that their continuum spectra are nearly indistinguishable. Therefore, 
the spectrum alone is insufficient to solve for idisk- As we now show, polarization can very 
effectively break the degeracy. 

5.3. Resolving the Spin-Inclination Degeneracy via Polcirization 

Measurement s 

Figure 5 shows the degree of polarization (P) as a function of photon energy for the 
degenerate models defined in Table 1. We see that (P) depends sensitively on the disk 
inclination. For the six models in Table 1, the disk inclination angle increases monotonically 
from model Ml to model M6, and we see that the degree of polarization of the disk radiation 
also increases monotonically. The difference in (P) between the extreme models is quite 
large. 

Figure 6 shows corresponding results for the polarization angle (ip). Once again wc see 
a large variation as a function of disk inclination. The variations are especially dramatic at 
photon energies on the order of several keV. Since the thermal emission from the disk peaks 
in this region of the spectrum (and photoelectric absorption is unlikely to be important), 
the effect could be easily measured (see § 6) . 

These results demonstrate that the degeneracy between the black hole spin and the disk 
inclination in the continuum blackbody spectrum of the accretion disk can be resolved by 
polarization measurements. Note that the degree of polarization, in particular, shows large 
variations with disk inchnation. To illustrate this fact, we show in Fig. 7 the run of (P) 
with photon energy for an accretion disk with M — 2 x 10^® gs~^ (similar to model M4 in 
Table 1) and a range of inclination angles idisk around a black hole of mass M = IOMq and 
spin parameter a* = 0.75. Over the range 30° < idisk ^ 80°, we see that (P) varies by about 
an order of magnitude. Thus, an observation of the polarization would provide an accurate 
determination of the inclination of the inner disk. A fit of the continuum spectrum would 
then give an accurate measurement of the spin of the black hole. 

As a by-product, such work would resolve the question of whether or not the inner disk 
in a black hole X-ray binary is aligned with the orbital plane (see § 2.2). Knowing the answer 
to this question would provide valuable constraints on models of core collapse and stellar 
black hole formation. 



-16- 



6. The Prospects for an Observational Test 

Is the inner accretion disk around stellar black holes warped away from the orbital plane 
(§2.2)? If so, is the warping severe enough to compromise the estimates of spin obtained using 
the continuum-fitting method (§§2.1 and 5.2)? In this section, we show schematically that 
we can expect definitive answers to these question in the near future through polarization 
measurements in the X-ray band. 

X-ray polarimetry is presently a hot topic, as evidenced by the 2004 X-ray Polarimetry 
Workshop at Stanford University^ and the upcoming conference to be held in Rome: The 
Coming Age of X-ray Polarimetry.^ Furthermore, several polarimetric missions have been 
proposed recently. We now briefly consider the performance of two very different instrument 
concepts developed within the severe constraints of a NASA SMEX-class payload. These 
are modest instruments because the mission costs for a SMEX (Small Explorer) are capped 
at $105 million (excluding the launch vehicle). 

The instrument considered by a team at the Goddard Space Flight Center (GSFC) 
is a broadband (1-10 keV) polarimeter, which measures the angular distribution of the 
tracks of the photoelectrons (Costa et al. 2001; Swank et al. 2004). The other instrument 
concept was developed by a team at the Smithsonian Astrophysical Observatory (SAO). 
It is a narrowband, Bragg-crystal instrument, which operates at 2.6 keV and 5.2 keV (E. 
Silver 2008, private communication) and is a modern version of the crystal polarimeter flown 
aboard OSO-8 (Weisskopf et al. 1978). The sensitivity of these two instruments is nearly 
identical and is summarized in Table 2. As indicated in the table, either instrument can, for 
example, detect polarization in a 1 Crab source at the 0.3% level in 1 day and at the 0.1% 
level in 10 days. 

Our focus here is the thermal dominant (TD) state of black hole binaries (§2.1), or an 
intermediate or steep power-law (SPL) state with a strong thermal component (McChntock 
& Remillard 2006), which we here refer to as a near-TD state. The dominance in this state 
of the thermal component makes it ideally suited for the determination of the inclination 
of the inner disk via polarimetry. Fortunately, nearly all of the many black-hole transient 
sources display a thermal dominant spectrum for weeks at a time during their outburst 
cycle. Furthermore, the sources are bright during this thermal phase, which makes them 
feasible polarimetric targets. Table 3 lists a selected sample of the more than 40 known 
transient and persistent black hole systems (McClintock & Remillard 2006; Remillard & 



^http://heasarc. gsfc.nasa.gov/docs/hcasarc/polar/polar. html 
''http://projects.iasf-ronia.inaf.it/xraypol/ 
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McClintock 2006; Orosz et al. 2007; Silverman & Filippenko 2008). The average intensities 
and outburst durations are restricted to those times when the source was in the TD or near- 
TD state. These data are only crude characterizations of the actual intensities/durations; 
see the references for precise information. 

As the data in Table 3 indicate, many sources have maintained a TD spectrum with an 
intensity of ~ 1 Crab for several tens of days. Therefore, a modest instrument like those 
described above would have the capability to measure the predicted polarizations of up to 
~ 5% (Fig. 7) with a sensitivity of ~ 0.1% for several transient (and persistent) sources 
during the course of a mission. Furthermore, polarization measurements of this precision 
could in many cases be repeated several times during the ~ 50-100 d thermal decay phase 
of a typical transient source. 

Even a single ~ 1-10 d observation would allow the inclination of the inner disk, for a 
wide range of inclinations, to be determined to roughly one degree (E. Silver 2008, private 
communication) . Meanwhile the inclination angle iorb of the orbital plane is already known 
for most of the sources in Table 3 to a few degrees (Charles & Coe 2006; Orosz et al. 2003). 
Thus, a simple comparison of the values of i^isk and iorh for a source in question will provide 
a stringent test of whether its disk is warped. 

7. Discussion and Conclusions 

In order to secure the measurement of black hole spin using the continuum-fitting 
method, it is essential to obtain an independent determination of the inclination of the 
inner accretion disk i^isk- We have shown that this appears to be entirely feasible via X-ray 
polarimetry. Our models predict polarizations of up to ~ 5% that vary monotonically with 
inclination, while sensitivities of ~ 0.1% arc achievable with quite modest SMEX-class po- 
larimeters. Such instruments are expected to be capable of routinely making measurements 
of ^disk over a wide range of inclinations with a precision of ~ 1°. Meanwhile, with full atten- 
tion it will be possible to achieve measurements of similar quality for the orbital inclination 
angle iorb, given recent advances in optical/NlR instrumentation (e.g., Vernet et al. 2007) 
and adaptive optics (e.g., van Dam et al. 2006). Once reliable values of these two inchnation 
angles have been obtained for a sample of black hole binaries, the question will simply be: 

is «disk ~ ^orb) is it not? 

There are of course significant hurdles that must be cleared. For example, the model of 
the disk atmosphere presented here includes only an approximate treatment of the effects of 
absorption. In the manner of Davis & Hubeny (2006), one must model the non-LTE effects. 
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Compton scattering, and the opacities due to ions of the abundant elements. Magnetohydro- 
dynamic models of thin disks in general relativity, which include radiation and polarization, 
are the ultimate theoretical goal, and progress is now being made on this front (Beckwith et 
al. 2008; Noble et al. 2008; Reynolds & Fabian 2008; Shafee et al. 2008b). Then, there are 
observational complications to consider. For example, even in the thermal dominant state 
there exists some remnant scattered coronal emission, which must be modeled. Some other 
examples of possible nettlesome sources of polarization include self-irradiation and reflec- 
tion in the disk, scattering by interstellar grains, and the effects of a global component of 
magnetic field in the disk. 

Clearly, the greatest hurdle is getting an X-ray polarimeter into space. Thirty years has 
passed since the tiny OSO-8 crystal instrument successfully measured the polarization of 
the Crab Nebula at 2.6 keV to be 19.2 ± 1.0% at a position angle of 156?4 ± 1?4 (Weisskopf 
et al. 1978). It is surely now time to open the polarimetric channel for serious exploration 
and discovery. In this paper we have shown how a polarimeter can secure the measurement 
of black hole spin. This is just one exciting application of X-ray polarimetry and there are 
numerous others. 

In the study of black holes, polarimetry furthermore promises to define the geometry 
of the emitting elements and constrain source models in decisive ways (§1). This is true, 
for instance, for modeling the rclativistic iron line (Miller 2007), or relativistic jets (Mirabel 
& Rodriguez 1999), or the mysterious power-law component that extends unbroken to McV 
energies (Grove et al. 1998), or the global oscillations that sometimes modulate up to one 
quarter of the total accretion power (McClintock & Remillard 2006). In modeling these 
and other black-hole phenomena, a central question is the geometry of the Comptonizing 
coronal source, which is vaguely and variously described as a sphere or a slab or a lamp post. 
Polarimetry will provide the best, and often only, clue to the actual geometry of the corona 
and other key structures (Meszaros et al. 1988; Blandford et al. 2002). 

This work was supported in part by NASA grants NNH07ZDA001N and NNX08AJ55G 
and NSF grant AST-0805832. We thank Eric Silver and Bill Forman for providing informa- 
tion on their instrument concept and Martin Elvis for a discussion on AGN. 

A. Relativistic Thin Accretion Disk Model 

We make use of a simple a-viscosity model to describe the disk (Shakura & Sunyaev 
1973), modified to include relativistic effects (Novikov & Thorne 1973; Page & Thorne 1974). 
Below, we first discuss the relativistic terms and then proceed to describe the disk model. 
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A.l. The Disk Height 

By equation (A8) of McClintock et al. (2006), the vertical gravitational acceleration in 
the comoving frame of the fluid, for small excursions z from the disk mid-plane, is given by 

9z = i^\j,z , (Al) 

where 



is the Newtonian angular velocity of the disk, R is the cylindrical radius, and 

e ^ i (1 - 4a.x-3 + 'ialx-') , ^ ^ ' ^ ^ " ^^3) 

The dimensionless function C, and the functions V, and Q used later below, are defined 
in Novikov & Thorne (1973) and Page & Thorne (1974). 

Hydrostatic equilibrium in the vertical direction gives 

I ^ g,{z = H) , (A4) 

where H is the height (vertical thickness) of the disk and Cg is the thermal sound speed of 
the gas at the disk mid- plane (both in the comoving frame). From equations (Al) and (A4), 
we obtain 

where the vertical disk frequency 

Qv = ^Ko^'^ . (A6) 



A. 2. The Flux Density and the Shear Stress of the Disk 

In this subsection we set G = c = 1 unless either of them or both appear explicitly in 
an expression. With these geometrized units, we have Rg — M and x — (R/My/"^. 

Page & Thome's functions / and Q (eqs. 15n and 35 of Page & Thorne 1974) are 
related by 

3 Q 
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Then, by equation (lib) of Page & Thorne (1974) and equation (Dll) of Li et al. (2005), 
we have the disk flux density (c.f. eq. 5.6.14b of Novikov & Thorne 1973) 



3GMM 1 / „ Cll'^x^ 



where ^2 is the angular velocity of the disk, E"^ and L'^ arc, respectively, the specific energy 
and the specific angular momentum of disk particles, g^^i = riti^M/VL^^ > is the torque 
applied at the inner boundary of the disk (at the ISCO), and ein = 1 — is the specific 
gravitational binding energy at the inner boundary (hereafter the subscript 'in' indicates 
evaluation at the disk inner boundary). 

By equation (11c) of Page & Thorne (1974), the vertically integrated coordinate-frame 
component of the shear stress tensor of the disk is 



■Mf + ^- p^Qi, 




27ri? ^ , 

J\> / /\/r\ 1/2 / /^l/2^3 \ 

(A9) 

where (p is the azimuthal coordinate in the disk. The vertically integrated comoving-frame 
component of the shear stress, W = 2Ht^j^ (hats indicate coordinate components in a frame 
comoving with the disk), is related to by (the last equation in page 428 of Novikov & 
Thorne 1973) 

= RBC'^/^V W . (AlO) 

Hence, equation (A9) agrees with equation (5.6.14a) of Novikov & Thorne (1973) in the limit 
of 77 = (zero torque at the inner boundary). Then we have 



Following Shakura & Sunyaev (1973) we describe the viscosity in the disk via a dimen- 
sionless parameter a and assume that t^f^ — (3/2)q;Pc) where is the pressure at the disk 
mid-plane, and the extra factor 3/2 is in order to be similar to the Newtonian equation. 
Then we obtain 
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Using 2HaPc — 2Hapcl — aECg and equation (A2), where E is the surface mass density of 
the disk, we then have 

ac^ ■ / (-.1/2 3 \ 

Note that the Newtonian hmit (a; 3> 1) of equation (A13) is 

Although in this paper we only consider the case 77 = 0, we keep 77 in the above formulae 
for generality. 



1/2 



(AM) 



A.3. One Zone Disk Model 

Wc make the standard one-zone approximation to describe the vertical structure of the 
disk (Frank et al. 2002). Thus, in terms of the mid-plane pressure Pc and density pc, the 
thermal sound speed and the surface mass density are given by 

cl = Pc/Pc, (A15) 
E = 2Hp^ . (A16) 

To allow for the effect of scattering and absorption opacity, we use the following ap- 
proximate formula to relate the escaping flux density F to the mid-plane temperature Tc, 

F = , (A17) 

(3/4)[r/2 + l/v^+l/(3r,b)] ' ^ ' 

where a is the Stefan-Boltzmann constant. Tab is the absorptive optical depth from the mid- 
plane to the surface of the disk, and r is the total optical depth. This formula is taken from 
Popham & Narayan (1995) and is based on the work of Hubeny (1990). Correspondingly, 
we use the following equation of state of the fluid. 



pckTc , 4(tT^ 



prriyt 3c 



r/2 + l/v^ 



r/2 + l/x/3 + l/(3rab) 



(A18) 



where the two terms correspond to gas and radiation pressure. Here, k is the Boltzmann 
constant, is the atomic mass unit, and fi is the mean molecular weight. We assume a 
fully ionized gas of solar composition, p = 0.603, and we approximate r and Tab as 

T = Tes + Tab , = 0.346 (E/2) , T^b = 6.6 X 10^^ p^T''/^ (E/2) , (A19) 
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where the coefficients are in cgs units. 

The disk solution is easily obtained using the above equations. At a given radius R and 
for an assumed value of E, equation (A13) allows us to calculate Cg, equation (A5) gives H, 
equation (A16) gives pc, equation (A15) gives Pc, equations (A18) and (A19) allow us to 
solve for Tc, and finally equation (A17) gives the radiation fiux F. Equating this estimate 
of the fiux to the correct value given in (A8) gives an algebraic equation with only one 
unknown, E, which can be solved numerically. Once we have E, we immediately obtain the 
values of r^s and Tab and this allows us to calculate the polarization suppression factor in 
equation (5). 



B. Propagation of the Polcirization Vector in Kerr Spacetime 

In this section we use geometrized units: G — c— 1. 

In terms of the Newman-Penrose orthonormal (ON) tetrad {P, n", m", m"} [where 
and are real vectors, m" and m" are complex, and = (m")* is the complex conjugate 
of m"], the Kerr metric tensor can be written as (Chandrasekhar 1983; Wald 1984)^ 

9ab = 2 [-l(anb) + m^ajrib)] , (Bl) 

where the use of parentheses '()' in subscripts denotes symmetrization of a tensor. The ON 
tetrad satisfies the orthogonality conditions 

lm = lm = nm = nm = 0, (B2) 

the normalization conditions 

l-n = — 1, m-m=l, (B3) 

as well as the null conditions 

ll = nn = mm = mm = 0, (B4) 
where 1 • m = Z"ma, etc. Then, for any vector we have 

f • f = -2[(f • l)(f • n) - (f • m)(f • m)] . (B5) 



*We adopt the signature convention of Misner et al. (1973) and Wald (1984), which is opposite to that 
of Chandrasekhar (1983). 



-23- 



Two Killing vectors and one Killing tensor in a Kerr spacetime lead to three conserved 
quantities along geodesies of particles and photons: the energy, the angular momentum 
about the symmetric axis of the black hole, and the square of the 'total angular momentum' 
(Carter 1968; Bardeen et al. 1972; Wald 1984). The three conserved quantities enable one 
to solve the timelike and null geodesies in a Kerr spacetime (Chandrasekhar 1983; Li et al. 
2005). 

The Kerr metric - like any other type D (II-ll) vacuum spacetime - possesses a conformal 
Killing spinor (Walker & Penrose 1970; Wald 1984) which enables one to determine the 
parallel propagation of polarization vectors along null geodesies in a simple manner. This is 
possible because of the Walker- Penrose theorem: If k is a null geodesic, afiinely parametrized, 
and f is a vector orthogonal to k and parallelly propagated along it, then, in a type-D 
spacetime, the quantity 

i^wp = 2 [(k • 1) (f • n) - (k • m) (f ■ m)] ^"'Z' (B6) 

is conserved along the geodesic (Walker & Penrose 1970; Chandrasekhar 1983), i.e., 

rVai^wp = . (B7) 

In equation (B6), ^2 is the only nonvanishing Weyl scalar in a Kerr spacetime 

,Tr / 2 Q 2 2 q\ .oMcOSi? . ^ 2 2 a\ ^ /roN 
W2 = ~"p~|6 V ~ '^^^ ^) ~ p~|6 V ~ ^ / ~ ^3 ' ^ ' 

where we have adopted Boyer-Lindquist coordinates {t^r^'d^ip), and 

= pp* = + a^cos^'(? . (B9) 



Two corollaries can be derived from the Walker-Penrose theorem. 

Corollary 1: Kwp is invariant under the transformation f — > f + ak, where a is an 
arbitrary function. 

Corollary 2: In a Kerr spacetime, the Walker-Penrose theorem implies that 

I^wpI' = [Q + (L, - aE^f] (f • f ) , (BIO) 

where £'00 = —h is the conserved energy-at-infinity, = k^^ is the conserved angular 
momentum about the axis of the black hole, and Q — k^ + Ll cot^ t!) — a^E"^ cos^ (When 
a = 0, Q + is the square of the total angular momentum.) 

Proof of Corollary 1 is straightforward by applying the identity k • k = and equation 
(B5). 



Proof of Corollary 2 is as follows. Following the proof of the corollary 2 in §60 of 
Chandrasekhar (1983), we have 



li^wpl' = 2 l^aP"^' (k ■ m)(k ■ m)(f ■ f) . (Bll) 



Then, defining 



S = -^ + asm^kt = ^^-asiwdE^ , (B12) 
sm 17 sm v 

T = = sign(A;^) ^JQ + a^E^ cos^ ^ - cot^ ^ , (B13) 

one finds that 

(k.m)(k.m) = -i-(5^ + r^) . (B14) 

It can be checked that 

S'^ + T'^ ^ Q+{L^- aE^f = constant , (B15) 

and 

|^2|-2/3^M-2/3|^|2_ (Big) 

Substituting equations (B14)-(B16) into equation (Bll), the identity in equation (BIO) is 
proved. 

Let us define 

XwP = {-M)-^'^ (K, + iK^) , (B17) 

where Ki and K2 are real constants. Then, substituting equation (B8) into equation (B6), 
we obtain 

+ iK2 = 2p* [(k • l)(f • n) - (k • m)(f • m)] . (B18) 

In terms of vector components in the Boyer-Lindquist coordinate system, we have 

Ki + iK2 = , / r { (r^ + a?) {Kft - hfr) + a (KU - Kfr) 

r + la cosv >■ ^ ^ 

"^sin^ [^^/v' - - « sin^ 1? {ktU - hft)] | , (B19) 
where we have used the fact that 

kaF = . (B20) 
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Equation (B19) is very useful since kt and are conserved quantities, and k^ can also be 
expressed in terms of conserved quantities (^'oo L^, and Q). Therefore, kr can be calculated 
using the equation g^'^kak^ — 0. 

To solve for the propagation of photon polarization, we take f to be the unit polarization 
vector, i.e. f = A/ A where A is the wave amplitude vector (Misner et al. 1973). The vector 
f satisfies equation (B20), and 



which states that the polarization vector is parallelly propagated along rays. Although for 
arbitrarily polarized radiation f is complex, for a linearly polarized beam it can be chosen 
to be a real vector (Misner et al. 1973). Hence, we assume that f is real, with 



The four-vector f is only defined to within a multiple of k since k satisfies the geodesic 
equation. Under a transformation f ^ f = f + ak, it is obvious that equations (B20) and 
(B22) are preserved. However, for f to be a solution of the parallel propagation equation 
(B21), a must satisfy /c^VaO; = 0, i.e. a must be a constant along the null geodesic. 

Because of the Walker-Penrose theorem, we do not need to solve the propagation equa- 
tion (B21) of the polarization vector explicitly. Instead, equations (B19) and (B20) can be 
used to determine the polarization vector. By Corollary 2, equation (B22) can be derived 
from equation (B19). Since /" is parallel propagated along the geodesic of k"-, faf"" and /a/c" 
are preserved along the geodesic. Hence, if at a point on the light ray equations (B20) and 
(B22) are satisfied, they are satisfied everywhere on the light ray. However, equation (B20) 
cannot be derived from the Walker-Penrose theorem so it is an independent equation. 

Equation (B19) is equivalent to two real equations. Then, together with equation (B20), 
we have three equations for the four components of The vector /" is determined up to an 
addition of a multiple of the wave vector A;" (Corollary 1). However, this degree of freedom 
does not affect physical measurements, since electromagnetic waves are transverse waves and 
a multiple of A;" only changes the component of f"" along the light propagation direction and 
the local time direction (Misner et al. 1973, Exercise 22.12). Hence, an uncertainty arising 
from a multiple of k"" does not prohibit us from making physical interpretations. In fact, 
we can make use of this gauge freedom to simplify the calculation by choosing a convenient 
form of Therefore, equations (B19) and (B20) are sufficient for solving for the physical 
components of After evaluation of the integral constants Ki and K2, equations (B19) 
and (B20) can be used to determine the component of /" at any point on the hght ray. 



k'^Vaf = , 



(B21) 



far = 1 ■ 



(B22) 
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Since we have assumed that f • f = 1, by equation (BIO), (B15), and (B17) we have 

Kl + Kl = S^ + T^ = Q+ (L, - aE^f . (B23) 

B.l. Evaluation of Ki and K2 on the Disk Plane 

On the disk plane we have ■»? = 7r/2, |pp = and ^2 — —M/r^. Equation (B19) is 
then reduced to 

Ki + iK2 = (r^ + a^) Ai + + i (A3 + 0^4) , (B24) 

where 

A, = r-' {krft - hfr) , (B25) 

A2 = r-' {KU - Kfr) > (B26) 

As = {hU - KM , (B27) 

A4 = r-' ihft - hU) . (B28) 

To evaluate the functions A1...A4 in terms of the local parameters of photons as they 
leave the disk, e.g., the photon energy and the propagation direction as measured by an 
observer corotating with the disk, we will make use of two local frames: the locally nonro- 
tating frame {e", e", e^, e^} which has a zero angular momentum, and the local rest frame 

|^(t)' ^(r)' ^(1?)' which corotates with the disk fluid (Bardeen et al. 1972; Novikov & 
Thorne 1973; Li et al. 2005). The two frames are related by a Lorentz transformation 



, - r [e^,) - v^el)] , (B29) 

= r [-v^el,^ + e^^)] , (B30) 

= e^,), (B31) 

= e^, (B32) 



where is the azimuthal velocity of a disk particle relative to the locally nonrotating frame, 

— 1 /2 

and r = (1 — v^) is the corresponding Lorentz factor. 

The direction of the velocity of a photon as it crosses the disk plane is specified by the 
normalized four-wavevector of the photon, = k°'/k^*^ — A;"/£'ioc, where k"' is the four- 
wavevector of the photon, and Eioc = fc^*^ = — /caC^^-j is the energy (frequency) of the photon 
measured in the local rest frame. The components of n"" in the local rest frame of the disk 
are 

n^*) = 1 , n^^'^ = - cos ^ , n^^) = sin ^ cos , n^'^^ = sin ^ sin , (B33) 
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where {9, 0) are spherical coordinates in the local rest frame, with the polar angle 9 measured 
from the normal to the disk, and the azimuthal angle measured relative to the vector e"^,-) 
along the disk radial direction (see Fig. 8).^ 

We write the polarization vector at the disk surface as 

r = fe^e + f% = (rcos^^cos0-/*sin0)e^,) + /sin^e^^) 

+ {f cos ^ sin + /'^ cos 0) e^^^ . (B34) 

By equation (B33), the wave vector is given by 

fc" = ^^loc [e"^) + sin 9 cos (f) e"^) - cos 9 e"^) + sin 9 sin e"<^)] . (B35) 

It can be checked that kaf"" = is satisfied. The condition = 1 implies that 

W +(/*)' = 1. (B36) 

Although with equation (B34) the polarization vector is defined up to an addition of a 
multiple of k"-, Corollary 1 indicates that this does not affect the values of Ki and K2. 



The disk angular velocity fl is related to the linear circular velocity by 

.2 \ 1/2 



(B37) 



where = + a^r(r + 2M), xo = {r"^^/ ^o)^^"^ (the lapse function in the equatorial plane; 
A = — 2Mr + a^), and luq = 2Mar/Ao (the frame dragging angular velocity in the 
equatorial plane). The specific angular momentum and the specific energy of disk particles 
are respectively 



(B38) 



With the above equations and the definition of the locally nonrotating frame (Bardeen 
et al. 1972; Li et al. 2005), we obtain 



E\oc 

-£"100 
AV2 



f E^cos9cos(f)- f 
-f L^cos9cos(f) + f 



rn 



Ac 



1/2 



sin 6^ + si 



sm ( 



f — j sin ^ + L''^ sin > , 



(B39) 
(B40) 



^Note, by definition, points inward to the disk. 
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^4 = E\oc { f 

Let us define 

X = 
Y = 



L^sin6' + r(^^ sin( 

sm ( 



+ cos^cos0S> , (B41) 

+ f'^rn[-^] cos^cos^V .(B42) 



1/2 



rr V A 

^0 



'''' [(r^ + a^) - aLt] = + - aA^^^^ , 



which satisfy 



Then by equation (B24) we have 

i^i = ^loc-^ [(/'cos^cos0-/^sin0)X-/^sin^F] , 

= -Eioc-^ [/sin^r + (/sin0 + /^cos^cos0)X] . 
It can be checked that equation (B23) is satisfied. 



(B43) 
(B44) 

(B45) 

(B46) 
(B47) 



B.2. The Limit at Infinity 

As r — oo, we have ~ 1 + M/r, and uu ~ 2Ma/r^. Expansion of the locally 
nonrotating frame at infinity leads to 



e*" ~ (1 



1 — 



r J \dt 



r J \dr 



+ 



1 



r [d^ ^ [ry r \dd 



rsini? \dip 



d 



O 



1 



dt J \r'^ J r \d(p 

dy " 

dr 



d 



J rsin'j? \d(p 



Omitting all corrections at or above the order of M'^ /r^, the metric takes the form 
^a6 ~ - f 1 - 



(B48) 
(B49) 
(B50) 
(B51) 



dtadtb + I 1 + — ) dvadn + r'^d^ad^b + sin^ Mipadipb . (B52) 
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We write the polarization vector /" in terms of the ON tetrad, 

r = fe: + fel + fel , (B53) 

where /'', f^, and must be finite. We have used the gauge freedom of to choose /* = 
at infinity. Then we have, as r — > oo, 

/t = 0, + U^rf\ U^rsin^f. (B54) 

Similarly, we can write the photon wave vector /c" in terms of the ON tetrad. Expressed in 
kn (which must be finite), we have 

\ r J \ T J r rsm(^^ 

= -ktel + /c^e" , as r ^ oo . (B55) 
The condition kah°- — then implies that (note, kr > 0) 

kr — —kt — Eoo , as r ^ oo , (B56) 
and equation (B20) leads to 

where we have omitted terms of the order r~^. 
Hence we have 

a{krf^-k^fr) ^ arsm.'dkrf + , (B59) 

k^U-KU ~ r(sm^k^f -k^f^+o{^^ , (B60) 

asm^^{ktU-k^ft) ^ arsm^^ktf + o(J^^ . (B61) 
Substituting these into equation (B19), we obtain in the limit r — > oo, 

K, + iK2 = ^(k^f + ^k^fA+asm'&krf 
rvj- \ Sin U J 

+^ [k^f - - a sin ^ k^A . (B62) 

\ sm v / 
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Substituting kr — —kt (eq. B56) into equation (B62), we have 

Ki + iK2 = - (rf + Sf^ - i (Sf - Tf^ , (B63) 

where S and T are defined by equations (B12) and (B13), respectively. Prom equation (B63) 
we can solve for and 

. ^ K,S + K,T ^ K,S - K,T 

where + is given by equation (B15) and is a constant. From equations (B23) and 
(B64), we then have 

/ -\ 2 n TV'S _|_ /^-s 

Since as r — > oo we have /* = and — (eq. B57), equation (B22) is satisfied. 



B.3. Solution of the Polarization Vector at Infinity 



ft = — jy^-^^^^^ {/[5(Xsin0 + Fsin^) -TXcos^cos. 



= — ' [f [-T{X sin + y sin 9) - SX cos 9 cos ( 



By equations (B46), (B47), and (B64) we have 
gAf{S' + T^) 

+f't'[T{Xsm(f) + Ysm9)+SXcos9cos(l)]} , (B66) 
gAf{S' + T^) 

+f'f'[S{Xsm(f) + Ysm9) -TXcos9cos(f)]} , (B67) 

where g = Eoo/Eioc is the photon redshift factor. Note that, on the right hand sides of the 
above equations, S and T are evaluated at the remote observer (hence — idisk, the disk 
inclination angle), but all other quantites are evaluated at the disk surface. 

The solutions given in equations (B66) and (B67) can be written in a matrix form 

cos $GR sin $GR 

-sin$GR COS$GR 





(B68) 



where 



cos$GR = = — 1/2°°^^ [5(Xsin0 + Fsin^) -TXcos^cos0] , (B69) 

gA^ {S'^ + T^) 

E r 

sin$GR = = — [T{Xsm(l) + Ysm9) + SX cos 9 cos (j)] . (B70) 

9^0 {S^ + 
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It can be checked that cos^ $gr + sin^ $gr = 1 is fulfilled. 
Define 

= cosV'oo , = sinV'oo , (B71) 

and 

f = cos -06111 , U = sinV'em ■ (B72) 

Then equation (B68) leads to 

cos ipoo = COS {ipeui - $gr) , sin ipoo = sin (V^em - $gr) , (B73) 
whose solution is 

ipoo = ^em - $GR + 2n7r , (B74) 

where n = 0, ±1,±2, .... Hence, the $gr defined by equations (B69) and (B70) represents 
the rotation of the polarization vector induced by the Kerr geometry and disk rotation. 

We define $gr to be the primitive rotation angle, which satisfies < $gr < 27r. Then 
by equations (B69) and (B70) we have 

^ f arccosfa; , if f „ > /t^„,-n 

= J ' ^ ~ . (B75) 

[ 27r — arccos , if ^j/ < 

Wc can check the nonrelativistic limit of this result in fiat spacctimc. For this purpose, let 
us take a = 0, f2 = 0, = 0, r = 1, L''" = 0, i?''' = 1, X = 1, (? = 1, and Aq = r^. Then 
we have X = and Y = 0. Since the light ray is not bent, wc have 6 = '&ohs = "^disk- 
So, we have = E'oor sin ^ sin 0, 5" = E'oorsin^, and T — sign (fc^^j^J \/ Q — cot^ 9 — 
—Eoorcos9cos(j). Then — 1, — 0, and hence the polarization vector is not rotated. 

The calculation of ^'gr and ipoo is simplified by using dimensionless variables. We define 
dimensionless quantities by symbols with a tilde: 

X = Xr^ , Y = Yr\ S = SrE^ , T = frE^ , Aq = Ar^ . (B76) 

Then, by equations (B69) and (B70), we have 

S(xsm(l) + Y sin 6^ - fX cos 6 cos 0] , (B77) 



6 = 



= ^ [f fxsin0 + ysin^) +,SXcos^cos0l , (B78) 
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where 



and 



--sin^?obs, (B80) 



r sin -i^obs r 

f = sign(A;^„,J V<9 + 0? cos^ i9ohs - cot^ ^?„bs , (B81) 



where A = L^/Eqo, and Q = Q/E'^ are independent of £'00 (Li et al. 2005). 
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Table 1. Parameters corresponding to the models shown in Figs. 4-6 



Model 


a* 


^disk 


M 


Ml 


0.998 


47.0 


0.48 


M2 


0.900 


60.0 


1.00 


M3 


0.830 


65.0 


1.40 


M4 


0.750 


70.0 


2.00 


M5 


0.630 


75.0 


3.20 


M6 


0.450 


80.0 


5.80 



Note. — Explanation of symbols: a*, dimension- 
less spin parameter of the black hole; idisk, inclination 
angle of the disk in degrees; M, mass accretion rate 
in units of 10^^ g s~^. Other parameters: mass of 
the black hole is lOM©; distance to the black hole is 
10 kpc; spectral hardening factor is 1.6. The inner 
boundary of the disk is at the ISCO and is assumed 
to have a zero torque boundary condition. 
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Table 2. Minimum Detectable Polarization (%)^ 



Intensity (Crab) 


0.001 


0.01 


0.1 


1 


Tobs = 1 d 


10 


3 


1 


0.3 


Tobs = 10 d 


3 


1 


0.3 


0.1 



^Nominal sensitivities accurate to ~ 20% (Swank et al. 2004; E. 
Silver 2008, private communication). 



Tabic 3. Selected Black Hole Binaries 





Coordinate 


Common 


Year/ 


Intensity'' 


Outburst^ 


References" 




Name 


Name/Prefix'' 


No. of Outbursts 


(Crab'^ ) 


Duration (d) 




1 


0620-003 


(A) 


1975/1 


15 


100 


Matilsky et al. 1976 


2 


1124-684 


Nova Mus 91 


1991/1 


1 


100 


Ebisawa et al. 1994 


3 


1543-475 


{4U) 


1971/4 


1 


100 


RM06 


4 


1550-564 


(XTE J) 


1998/5 


2 


100 


RM06 


5 


1650-500 


(XTE J) 


2001/1 


0.5 


100 


MR06 


6 


1655-40 


(GRO J) 


1994/3 


1 


300 


RM06 


7 


1659-487 


GX 339-4 


1972/12 


0.5 


200 


RM06 


8 


1705-250 


Nova Oph 77 


1977/1 


1 


50 


Watson et al. 1978 


9 


1859+226 


(XTE J) 


1999/1 


0.5 


100 


RM06 


11 


2000+251 


(GS) 


1988/1 


3 


100 


Tsuncmi et al. 1989 


12 


1630-472 


(4U) 


1971/16 


0.2 


100 


MR06 


13 


1743-322 


(H) 


1977/4 


0.5 


100 


RM06 


14 


0538-641 


LMC X-3 


Persistent 


0.03 




MR06 


15 


0540-697 


LMC X-1 


Persistent 


0.02 




MR06 


10 


1915+105 


(GRS) 


Quasi-persistent* 


1 




McClintock ct al. 2006 


16 


1956+350 


Cyg X-1 


Persistent^ 


1 




Wilms et al. 2006 



"A prefix to a coordinate name is enclosed in parentheses. 

''Crude estimate of average intensity in TD or near-TD state; see references. 

= 1 Crab Si 2.8 X 10"* erg cm"^ s"^ (1-10 keV) for a Crab-like spectrum with F = 2.1. 

Crude estimate of time in TD or near-TD state; see references. 
"RM06 = Remillard & McClintock 2006; MR06 = McClintock & Remillard 2006. 
*In outburst since discovery in 1992; occasionally observed in TD state (McClintock et al. 2006). 
^Observed in a soft, near-TD state for extended periods (e.g., Wilms et al. 2006). 
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Fig. 1. — Degree of polarization as a function oi /j, — cos 9 for a semi- infinite scattering 
atmosphere, where 9 is the angle between the line-of-sight and the normal to the plane of 
the atmosphere. 
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Fig. 2. — Polarization angle and degree of polarization as a function of photon energy for 
disk models with various disk inclinations ijisk and black hole spins a*. Following Connors 
et al. (1980), we have set the black hole mass M = 9Mq and the mass accretion rate 
M = 7 X 10^'' g s~^. The spectral hardening factor is /coi = 1-6. The inner boundary of 
the disk is assumed to be at the ISCO and to have a vanishing torque. The various hues 
correspond to different models, as indicated. All the models computed by Connors et al. 
(1980) have been recomputed and are shown, along with a number of additional models. 
The dotted lines correspond to the Newtonian hmit. 
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Fig. 3. — Continuum disk spectra computed with the pubhcly-available relativistic accretion 
disk model KERRBB (Li et al. 2005) for the three disk models shown in Fig. 7 of Connors et 
al. (1980). The distance to the source is assumed to be 10 kpc and the spectral hardening 
factor is taken to be 1.6. 
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Fig. 4. — Disk continuum spectra computed with kerrbb for the models M1-M6 defined in 
Table 1 and identified by the disk inclination angle i. Although the black hole spin parameter 
is very different in the various models, the continuum spectra are nearly indinstinguishable. 
This shows that the X-ray continuum spectrum alone cannot be used to determine both the 
spin parameter and the disk inclination. 
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Fig. 5. — Degree of polarization versus photon energy for the set of degenerate disk models 
M1-M6 defined in Table 1 and identified by the disk inchnation angle i. The computations 
were done using an extended version of kerrbb. In contrast to Fig. 4, where the continuum 
spectra of these models were identical, here the models are seen to have widely different 
degrees of polarization. 
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Fig. 6. — Similar to Fig. 5, but showing the angle of polarization. Once again the models 
are easily distinguishable. 
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Fig. 7. — Degree of polarization versus photon energy for a series of accretion disk models 
with varying inclination angle idisk- All the models have M = 2 x 10^^ gs~^, M = IOMq and 
= 0.75. The large variation in the degree of polarization with i^isk nieans that it should 
be possible to determine the disk inclination accurately with polarization data. 
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Z 




X 



Fig. 8. — The coordinate system on the disk plane. The disk center is located at O. 
{e"^), e"^), e"^), e"^)} is a local rest frame of the disk, with not shown. A hght ray is 

— * 

emitted by a disk particle at P, propagates along a direction labeled by k. The plane PAB 
is a meridian plane at P, perpendicular to the disk plane. In the spherical coordinate system 
{0,(j)), 9 is measured from the disk normal — e^^) in the meridian plane, and is measured 
from e"^^ in the disk plane. The unit vectors (in the meridian plane and perpendicular to 

k) and (in the disk plane and perpendicular to PB) are shown. 



